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The dynamics of an exciton-polariton superfluid resonantly pumped in a semiconductor micro- 
cavity are investigated by mean-field theory. Modulational instability develops into crystalline order 
and then ordered and disordered states alternately form. A supersolid-like state is also found, in 
which superflow coexists with crystalline order at rest. 
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I. INTRODUCTION 

Transitions between ordered and disordered states can 
be controlled by changing the parameters of a system. 
For example, a liquid can be solidified and then melted 
back into a liquid by controlling the temperature or pres- 
sure. Another example is a flow pattern in the wake of an 
obstacle. A periodic vortex street can change into turbu- 
lence by increasing the Reynolds number and turbulence 
can transform back into a vortex street by decreasing the 
Reynolds number. In a large complex system, order and 
disorder can alternately appear in a partial system, such 
as clouds in the sky. However, in a simple laboratory 
system, such repeated transitions between ordered and 
disordered states usually do not occur unless an external 
parameter is changed 1 . Ordering phenomena or pattern 
formations occur in a wide variety of nonequilibrium dis- 
sipative systems 2 -, such as Benard convection cells 3 * and 
chemical patterns in reaction diffusion systems^. How- 
ever, once a pattern has been generated in these systems, 
it is preserved if the parameters are fixed, although it may 
vary with time. By contrast, we here propose a system in 
which crystalline order and disorder alternately emerge 
despite external parameters being fixed. 

The system that we consider is an exciton-polariton 
superfluid in a planar semiconductor microcavity^i 
There have been many recent experimental studies of 
this system, including: bistable properties^, Bose- 
Einstein condensatio n 10 ^ 1 , Bogoliubov excitations^, 
quantized vortice a 13 i 14 , superfluidity 15 , quantum hydro- 
dynamics 16,17 , and pattern formation 18,19 . Of the var- 
ious ways to create a polariton superfluid, we con- 
sider a system of polaritons pumped by a quasi-resonant 
continuous- wave laser. Since polaritons have a short life- 
time of ~ 10 ps, which is comparable to the time scale of 
the dynamics, the system is a nonequilibrium dissipative 
system with continuous feeding. It is thus suitable for 
studying nonlinear phenomena far from equilibrium. 

In this paper, we investigate pattern formation and 
destruction dynamics in an exciton-polariton superfluid. 
The polaritons are first pumped to the lower branch of 
the bistability 8,9 , which is an unstable stationary state 20 . 
Modulational instability then arises and a hexagonal 
crystalline pattern develops, which eventually collapses, 
destroying the order. However, the pattern reappears 



from the disordered distribution and grows again. Subse- 
quently, the ordered and disordered distributions emerge 
alternately. We also show that the crystalline order can 
remain at rest even in the presence of superflow. 

This paper is organized as follows. Section [TT] formu- 
lates the problem, Section IIIII demonstrates numerical 
results, and Section HV1 gives conclusions of this study. 



II. FORMULATION OF THE PROBLEM 

Excitons in quantum wells and photons confined be- 
tween Bragg mirrors are assumed to be described by 
macroscopic wave functions ipx and iJjq , obeying the two- 
dimensional nonlinear Schrodinger equations given by 20 
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where Hq and jj are respectively the free Hamiltonians 
and the decay rates of an exciton (j = X) and a pho- 
ton (j = C), is the Rabi frequency, and g is the 
exciton-exciton interaction coefficient. The free Hamil- 
tonians are approximated as = Tuuj$ and Hq = 
hujQ — h 2 V 2 1 \2mc) , where rac is the in-plane effective 
mass of a cavity photon. Diagonalizing the first and sec- 
ond terms on the right-hand side (rhs) of Eq. (pQ), we 
obtain the dispersion relations uu±(k) for the upper and 
lower free polaritons. The last term on the rhs of Eq. (fTa|) 
describes pumping by an external laser beam with a pro- 
file F(r), in-plane wave vector fe p , and frequency uj p . 

We numerically solve Eq. (pQ) using the pseudospectral 
method 21 . As an initial state, small white noise is set 
to iJjx and ipc to exclude exact numerical symmetry. In 
the following calculations, we only consider the case of 
normal incidence of the pump beam on the sample, i.e., 
k p = 0. The detuning of the pump frequency from the 
lower polariton is defined as S = u p — cj_(0). We use the 
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FIG. 1: (color online) (a)-(i) Snapshots of the density profile 
\ipc\ 2 of photons. The pump parameters in Eq. (j2) are hFo = 
8.97 meV and d — 50 /im. The color scale in each image is 
normalized by the peak density, where the unit of density is 
/im" 2 . The field of view of each panel is 100 x 100 /mi. (j)-(l) 
Fourier-transformed profiles | f i/jce~ lk r dr\ 2 , where the field 
of view is 4 x 4 /im -1 . (m) Time evolution of Ac = f \^c\ 2 dr 
(red solid line) and Nx = f \ipx\ 2 dr (green dashed line). The 
vertical dotted lines show t = 400, 550, and 610 ps. See 
supplemental material. 



Gaussian profile of the pump beam as 



F(r) = F exp [-(x 2 + y 2 )/d 2 ] 



(2) 



where Fo is proportional to the peak intensity and d is 
the 1/e width. The parameters are rac = 2 x 10 _5 m e 
with m e being the electron mass, ftf^R = 5 meV, g = 
0.01meV/im 2 , and hS = 1.53 meV. For simplicity, we 



III. NUMERICAL RESULTS 

Figures Q](a)-[2(i) show the dynamics of the density 
profile \ipc(r)\ 2 of photons (\ftx(r)\ 2 is always similar to 
\ftc(r)\ 2 ', it will not be shown henceforth). First, the po- 
lariton density grows in the region of r < d and reaches a 
steady state with a Gaussian profile (Fig. [2(a)), reflect- 
ing the pump profile in Eq. (j2j). Modulational instability 
then arises around the center of the pump region and 
a pattern develops (Figs. [2(b) and [2(c)). At this stage, 
the pattern has no crystalline symmetry. At t ~ 500 
ps, crystalline order emerges around the center and it 
grows outward to form a triangular lattice (Figs. [2(d) 
and [2(e)). The wave function has almost the same phase 
at each density peak. This ordered state is not a steady 
state and the density of each peak continues to increase. 
When the peak density reaches some critical value, the 
triangular lattice abruptly collapses, destroying the crys- 
talline order (Figs. [2(f) and [2(g)). The disturbance de- 
cays with time and a similar modulational pattern to that 
in Fig. [2(c) re-emerges (Fig. [2(h)), which is followed by 
the reappearance of crystalline order (Fig.[2(i)). The col- 
lapse and restoration of crystalline order are subsequently 
repeated. 

Figure GJni) shows the time evolution of the norms 
N c = J \ipc\ 2 dr and N x = J \ftx\ 2 dr. After the initial 
fluctuations for t < 100 ps, Nq and Nx reach plateaus 
corresponding to the steady state in Fig. [2(a). The 
Fourier-transformed profile of this steady state localizes 
within \k\ < 0.1 /im -1 (data not shown). As the modu- 
lational instability grows, a ring-shaped profile emerges 
at \k\ ~ 1.3 /im" 1 ~ 27r/(5/mi) (Fig. [Tp)), which re- 
flects the pattern in Fig. [2(c). The anisotropic shape of 
the ring depends on the small noise added to the initial 
state. As the crystalline order grows, Nq and Nx in- 
crease rapidly and the ring in Fourier space changes to 
hexagonal peaks (Fig.QJk)). Since the peaks in Fig.QJk) 
are closer to the origin than the ring in Fig. DID, the pe- 
riod of the triangular lattice in Fig. [2(e) is larger than 
the modulation wavelength in Fig. [2(c). The sudden re- 
duction in Ac and Ax corresponds to the collapse of the 
triangular lattice, where the six hexagonal peaks change 
to a disordered distribution in Fourier space (Fig. [2(1))- 
The disordered components are subsequently cleaned up 
and the ring shape grows, which is followed by hexagonal 
peaks. 

To qualitatively understand the behaviors shown in 
Fig. [2 we consider the lower-polar iton approximation. 
Neglecting the upper polariton branch, we can reduce 
Eq. Q to2i 
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where ft LP = X LP ^x + Clp^c, rn LF = 2ra c , #lp = 
|^Lp| 4 # 5 and Flp = C^pF are the macroscopic wave 
function, the effective mass, the interaction coefficient, 
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and the pump intensity of the lower polar itons, respec- 
tively. Xlp and Clp are the Hopfield coefficients. For 
simplicity, we consider the case of homogeneous pumping 
and the polariton wave function can be divided into 

^Lp(r,t) = e-^ t [0 ss +^(r ) t)]. (4) 

Substituting Eq. (jlj into Eq. (j3j), we find that the homo- 
geneous steady state (j) ss is determined by 

(ns - #lp|0ss| 2 + *fry/2) = fi^LP, (5) 

which will have multiple solutions of </> ss for a certain 
range of Flp if S/j > y/3/2. For the parameters in 
Fig. [TJ (5/7 23 and Eq. ([5]) has three solutions for 0.8 
meV < hF^p < 14 meV. The steady state in Fig. QJa) 
corresponds to the state with the lowest density of the 
three solutions. 

We use Bogoliubov analysis to investigate the stability 
of the steady state. The first order of 5(f) in Eq. (j3]) gives 
the Bogoliubov-de Gennes equation. For an excitation 
of the form 5<j>(r,t) = u k e~ iujt+ik - r + v ^e iujH - ik - r , the 
Bogoliubov-de Gennes equation can be diagonalized to 
yield the excitation frequency: 

1 i 

uj = -i- ± V(^fc - S + u 88 )(u k - S + 3cj ss ), (6) 

where uo k — ?ifc 2 /(2mLp) and co> ss = ^lp|0ss| 2 /^- When 
the imaginary part of uj is positive, the corresponding 
mode grows exponentially. From Eq. (j6]), the condition 
for this modulational instability to occur is found to be 
uj ss > 7/2 and S - 2uj ss + (uj 2 s - 7 2 /4) 1 / 2 > 0, which 
is satisfied when r < 20 /im in the situation shown in 
Fig. [T{a). For uj k = 5 — 2cj ss , Eq. (j6j) has the largest 
imaginary part, which corresponds to the most unstable 
mode. For the parameters in Fig. [U^a), the most unsta- 
ble wavelength is estimated to be ~ 7 /im, which is in 
reasonable agreement with the modulation in Fig. QJb). 

The hexagonal crystalline order shown in Figs. [TJe) 
and Q] (i) is a manifestation of the nonlinear effect. To 
see this, we simplify the problem using a wave function 
with the form 

Mr, t) = e-^Vss[l + c(t)/(r)], (7) 

where f(r) = ^cosfe^ • r. Substituting Eq. ([?)) into 
Eq. ((3]) and extracting terms proportional to /(r), we 
obtain 

= (cj fe -cj p - ^7/2) c + cj ss [2c + c* + A(2\c\ 2 + c 2 ) 
+5|c| 2 c], (8) 

where A and 5 are constants depending on f(r). The 
0(c) terms in Eq. (JHj), which give Eq. (|6j), are responsible 
for the exponential growth of c and we restrict the wave 
number \ki\ to be the most unstable one. If A ^ 0, 0(c 2 ) 
terms will become important as c grows. Equation (|8j) 
gives 

i.| c |2 = _ 7 | c |2 + i0Jss{c 2 _ c *2 } + i0JssA \c\\c - C*), (9) 
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FIG. 2: (color online) (a)-(c) Snapshots of the density profile 
|^c| 2 for d = 60 /xm. The insets show Fourier- transformed 
profiles. The field of view of the main panels is 100 x 100 /im 
and that of the insets is 4 x 4 /im -1 . (d) Spatiotemporal image 
of the density on the dashed line in (a), (e) Time evolution 
of the norms Nc (red solid line) and iVx (green dashed line) . 
In (a)-(e), all the parameters except d are the same as those 
in Fig. [T] See supplemental material 22 . 



which indicates that the solution satisfying iA(c—c*) > 
dominates the exponential growth, breaking the sym- 
metry with respect to c — c in the linear part of 
Eq. (jSJ). The simplest form of f(r) with nonzero A is 
f(r) — cos ki-r+cos k2-r+cos k3-r with /ci+&2+&3 = 0; 
i.e., the three wave vectors form a regular triangle. In 
this case, A = 1 and B = 15/4. Thus, the hexagonal or- 
der dominates the growth as the modulation amplitude 
grows. The term of B in Eq. (j8]) effectively shifts the res- 
onance by u ss B\c\ 2 . This is why the peak wave number 
decreases (Fig. Hfk)) as the hexagonal pattern grows. 

We return to the numerical simulations of Eq. (pp). 
When the peak density grows and exceeds some criti- 
cal density, the subsequent behaviors can be classified 
into three types. The first one is the behavior shown 
in Fig. [TJ where the hexagonal order and disorder alter- 
nately emerge. The second behavior is entrainment into 
the large density for a large pump intensity Fo. When 
the local peak density reaches the uppermost branch of 
Eq. (j5j) , the whole system is entrained into the uppermost 
branch and the dynamics settle down there. The third 
behavior occurs when the width d of the pump beam is 
large, which is shown in Fig. O After the hexagonal pat- 
tern is formed and destroyed at t ~ 400 ps in a manner 
similar to Figs. GIa)-[3Jf), the collapse and restoration of 
the peaks occur locally and sporadically, which results in 
the pattern shown in Figs. Efa)-[2fc). This is in contrast 
with the dynamics in Fig. [TJ where the order-disorder os- 
cillation occurs synchronously in the whole system. Fig- 
ure EJd) shows a spatiotemporal image of the density 
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FIG. 3: (color online) (a)-(c) Snapshots of the density \i^c\ 2 
and phase arg^c profiles for Laguerre-Gaussian pumping. 
The pump parameters in Eq. (|10p are TiFq — 21 meV and 
d — 50 /mi. The field of view of each panel is 120 x 120 
/im. (d) Time evolution of the normalized angular momenta 
Lj — — i f ipjdeipjdr/ f \ipj\ 2 dr of photons (j = C, red solid 
line) and excitons (j = X, green dashed line). The insets 
show Fourier-transformed profiles, where the field of view is 
4x4 jLtm -1 . Except for the pump function, the parameters 
are the same as those in Fig. [1] See supplemental material 22 . 

|^c| 2 on the dashed line in Fig. [2fa). Each density peak 
blinks periodically with a period ~ 4 ps. This may be 
related with the blinking phenomenon observed in ex- 
periments^. Since the peaks grow and collapse at various 
places, the norms Nq and TVx fluctuate around their av- 
erages, as shown in Fig. Efe). There is also a case in 
which several cycles of the synchronous order-disorder 
oscillations shown in Fig. [1] occur, followed by the local 
oscillations shown in Fig. [2j 

We next examine polariton pumping with an angular 
momentum using a Laguerre-Gaussian beam. The pump 
function is given by 

F(r) = Fo^-exp [-(x 2 + y 2 )/d 2 ] . (10) 

The pumped polaritons form a ring-shaped density pro- 
file with a phase = arg(x + iy) and the modulational in- 
stability grows along the ring, as shown in Fig.[3fa). The 
triangular lattice initially has some mismatches along the 
ring, as shown in Fig[3fb), since the lattice grows from 
different points on the ring. Subsequently, in a man- 



ner similar to Fig. [2j the peaks collapse and revive lo- 
cally (Fig. EJc)). The angular momenta and the Fourier- 
transformed profiles shown in Fig. [3^d) exhibit an inter- 
esting behavior. Although the angular momentum per 
particle is ~ Ti on average, from which the angular ve- 
locity for a typical radius R ~ 40 fim is estimated to be 
~ h/ '(rriLpR 2 ) ~ 1.8 rad/ns, the hexagonal pattern in 
Fig. [3jd) does not rotate; rather, it drifts slightly in the 
opposite direction (clockwise). This indicates that the 
hexagonal crystalline order is almost irrotational despite 
the rotation of the superflow^; this may be regarded as 
a supersolid-like behavior. 

In experiments, a disorder potential in a sample typ- 
ically has an energy of at least ~ 0.1 meV and a 
wavelength of ~ 10 /im^. We have numerically con- 
firmed that dynamics similar to Figs. [THU are observed 
even when the disorder potential is added to the rhs of 
Eq. ([Taj). Under continuous pumping of polaritons such 
as in the present situation, it is experimentally difficult to 
perform time-resolved measurements. Time-integrated 
imaging will reduce the hexagonal order, since its direc- 
tion is variable (see Figs. QJe) andQJi)) because of the 
rotational symmetry of the Hamiltonian. However, in re- 
ality, the disorder potential breaks the rotational symme- 
try and fixes the direction of the hexagonal order, which 
enables us to use time-integrated imaging to identify the 
hexagonal order. The order-disorder oscillation shown in 
Fig. [T] may also be detected by taking spatial and tem- 
poral correlations^. 



IV. CONCLUSIONS 

We have investigated the dynamics of an exciton- 
polariton superfluid coherently pumped in a semiconduc- 
tor microcavity. Under appropriate pumping, a hexago- 
nal crystalline order is formed, destroyed, and restored; 
this cycle is repeated many times (Fig. [I]). Such a sponta- 
neous oscillation between ordered and disordered states 
is a unique phenomenon in a simple nonequilibrium open 
system. We have also shown that superflow can coexist 
with crystalline order at rest (Fig. [3j), which is reminis- 
cent of supersolidity. The continuously pumped exciton- 
polariton system thus affords a quantum dissipative sys- 
tem that is promising for studying novel dissipative struc- 
tures. 
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